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LoCuSS: the connection between brightest cluster galaxy activity, 
gas cooling and dynamical disturbance of X-ray cluster cores 
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ABSTRACT 

We study the distribution of projected offsets between the cluster X-ray centroid and the 
brightest cluster galaxy (BCG) for 65 X-ray selected clusters from the Local Cluster Sub- 
structure Survey (LoCuSS), with a median redshift of z = 0.23. We find a clear correlation be- 
tween X-ray/BCG projected offset and the logarithmic slope of the cluster gas density profile 
at 0.04r5oo (a), implying that more dynamically disturbed clusters have weaker cool cores. 
Furthermore, there is a close correspondence between the activity of the BCG, in terms of 
detected Ha and radio emission, and the X-ray/BCG offset, with the line emitting galaxies all 
residing in clusters with X-ray/BCG offsets of ^15kpc. Of the BCGs with a < —0.85 and 
an offset < 0.02r5oo, 96 per cent (23/24) have optical emission and 88 per cent (21/24) are 
radio active, while none has optical emission outside these criteria. We also study the cluster 
gas fraction (/gas) within r^Qo and find a significant correlation with X-ray/BCG projected 
offset. The mean /gas of the 'small offset' clusters (< 0.02r5oo) is 0.106 ± 0.005 (cr = 0.03) 
compared to 0.145 ± 0.009 (cr = 0.04) for those with an offset > 0.02r5oo, indicating that the 
total mass may be systematically underestimated in clusters with larger X-ray/BCG offsets. 
Our results imply a link between cool core strength and cluster dynamical state consistent 
with the view that cluster mergers can significantly perturb cool cores, and set new constraints 
on models of the evolution of the intracluster medium. 

Key words: galaxies: clusters: general - X-rays: galaxies clusters - cooling flows - galaxies: 
evolution - galaxies: elliptical and lenticular, cD 



1 INTRODUCTION 

The first ranked, or brightest cluster galaxies (BCGs) are amongst 
the brightest of all galaxies, accounting for around 5-10 per cent of 
the total light in massive clusters (Lin & Mohr 2004). Their prop- 
erties are closely related to those of the host cluster (Edge 1991; 
Lin & Mohr 2004) and they typically lie at the bottom of the poten- 
tial well. In relaxed clusters, this location is frequently positioned 
within a 'cool core' in the intracluster medium (ICM), where gas 
is capable of condensing out of the hot phase, to form stars (e.g. 
Crawford et al. 1999; Rafferty, McNamara & Nulsen 2008) or fuel 
accretion onto a central supermassive black hole (e.g. Best et al. 
2007). In this configuration the BCG lies at a interface which is 
crucial to understanding the role of feedback in galaxy and cluster 
evolution. 

The unique character of BCGs allows them to be used as 
important diagnostics of the host cluster. In particular, the pro- 
jected offset between the X-ray peak and the BCG is sensitive 
to the cluster dynamical state (Katayama et al. 2003), which is 
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analagous to the offset between the peaks of the gravitational lens- 
ing mass map and the X-ray emission as an indicator of distur- 
bance (Smith et al. 2005). Furthermore, the presence of Ha emis- 
sion from ongoing star formation (e.g. Heckman 1981 ; Peres et al. 
1998; Crawford et al. 1999; Edwards et al. 2007; Cavagnolo et al. 
2008) and radio emission from active galactic nuclei (AGN) activ- 
ity (e.g. Peres et al. 1998; Best et al. 2007; Cavagnolo et al. 2008; 
Mittal et al. 2008) probe the thermodynamic state of the hot gas 
in the cluster core. For example, it has recently been discovered 
that star formation in BCGs only occurs when the entropy of the 
ICM falls below a critical threshold (Voit et al. 2008; Rafferty et al. 
2008). Moreover, almost all clusters with star-forming BCGs lie 
above the X-ray luminosity-temperature relation (Bildfell et al. 
2008), demonstrating the close link between the BCG and global 
cluster properties. 

With the advent of Chandra it is possible to probe cluster 
gas properties on the scale of BCGs out to intermediate redshifts 
and with its extensive archive of observations, this type of analy- 
sis can be applied to large numbers of clusters. In this paper we 
use Chatidra data to explore the connection between cooling in 
cluster cores and the location and activity of the BCG for an X- 
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ray selected sample of 65 clusters drawn from the Local Cluster 
Substructure Survey^ (LoCuSS). This is a morphologically unbi- 
ased sample of ~ 100 X-ray luminous galaxy clusters selected from 
the (e)BCS and REFLEX ROSAT All Sky Survey (RASS) cata- 
logues (e.g. see Zhang et al. 2008; Okabe et al. 2009), containing 
all clusters down to the RASS flux limit, within the LoCuSS red- 
shift, declination and Galactic column cuts. Throughout this paper 
we assume Ho = 70 km s"^ Mpc^^ , fim ~ 0.3 and Q,a = 0.7. 
Errors are quoted at a 68 per cent confidence level. 



2 SAMPLE SELECTION AND DATA ANALYSIS 

We select all clusters with available Chandra data that satisfy the 
LoCuSS selection function of 0.15 < z < 0.3, -70° < (5 < 70°, 
riH < 7 X 10^" cm^^. This produces a total of 66 clusters, from 
which we discard ZwCl 1 309. 1 +22 1 6 because the data on this clus- 
ter are too shallow for our purposes. The median redshift is z = 0.23 
with a corresponding scale factor of 3.7 kpc per arcsecond. Eigh- 
teen of the 65 cluster observations are drawn from our own Chan- 
dra programs in Cycles 9 and 10 (PIDs: 09800732, 10800565), the 
remaining observations are from the archive. 

2.1 X-ray data analysis 

Chandra data were reduced and analysed according to the 
procedure described in Sanderson, O'SuUivan c& Ponman (2009). 
Briefly, the Chandra data were reprocessed using CIAO version 3.4 
and incorporating CALDB version 3.4.2, to produce flare cleaned 
and point source removed level 2 events files. Corresponding blank 
sky background datasets were also produced and matched in nor- 
malization to the cluster events in each case, to account for vari- 
ations in the particle-dominated high energy background. No ad- 
justment was made to the background to allow for any variation 
in soft Galactic foreground emission compared to each cluster ob- 
servation. However, following Sanderson, Ponman & O'SuUivan 
(2006), the galactic absorbing column was fitted as a free parameter 
in the spectral modelling to allow for any differences in inferred low 
energy absorption associated with soft emission excesses or cali- 
bration uncertainties. For some clusters, it was necessary to freeze 
the absorbing column at the galactic HI value, since unfeasibly low 
values were otherwise obtained in many of the annular bins. 

A series of annular spectra were extracted for each cluster, 
comprising between 6 and 25 radial bins in total, depending on 
data quality. These annular spectra were then fitted with an ab- 
sorbed APEC model using the PROJCT scheme in XSPEC version 
1 1.3.2 to yield deprojected gas temperature and density profiles for 
each cluster. Weighted response matrix files were used for each 
spectrum, and the fitting was performed in the energy range 0.5- 
7.0 keV. In order to provide more reliable error estimates on the 
annular spectral measurements, we performed 200 Poisson real- 
izations of the best-fit PROJCT model, using the same background 
and response files (Sanderson, in preparation). Each realization was 
treated like the original data and fitted with the PROJCT scheme to 
produce a suite of simulated measurements in each bin. The errors 
in each bin were then calculated from the median absolute devi- 
ation of these simulated values, as a robust estimator of the stan- 
dard deviation that is well-suited to heavy-tailed distributions (see 
Beers, Flynn & Gebhardt 1990, for example). 

^ http://www.sr.bham.ac.uk/locuss 



We note that a potential issue with our Chandra analysis is 
the possibility of a modest bias in temperature estimates for hotter 
clusters (>4-5 keV) resulting from errors in the Chandra response 
matrix (as described in Sun et al. 2009, and references therein), 
which could lead to overestimates of the temperature. To assess 
the magnitude of the effect this could have on our analysis, we 
have reanalysed one of the cluster datasets using CIAO version 4.1, 
incorporating CALDB version 4.1.2. We have selected Abell 1835 
(obsid 6880), since this is a deep observation (120ks) of a bright, 
massive cluster, where the impact of calibration changes is likely 
to be greatest. Using the new calibration data, we find that rsoo 
changes to 1432 ± 60 kpc (from 1506 ± 57, in our original analy- 
sis), with a corresponding gas fraction within rsoo of 0.128±0.011 
(0.110 ± 0.008 originally). As expected, the newer CALDB results 
in a lower 7-500 and hence higher gas fraction (due to both a sys- 
tematic lowering of the temperature leading to a lower total mass, 
as well as an increase in the gas mass measured from the emissiv- 
ity), but the effect is only ~5 per cent for rsoo and 16 per cent for 
the gas fraction; corresponding to 1.2 and 1.6(j, respectively. Since 
only two of our clusters are cooler than 4keV, such a systematic 
shift will be similar across the whole sample and so is not likely to 
impact our results significantly. 

The centroid for the spectral profile analysis was used to de- 
termine the projected offset between the centre of the cluster X-ray 
halo and the BCG. We determined this position in a similar fashion 
to Maughan et al. (2008), iteratively refining the cluster X-ray sur- 
face brightness centroid within an aperture of 2-3 arcmin initially, 
then repeating the centroiding within an aperture of ~1 arcmin. In 
addition, we also determined the coordinates of the X-ray peak as 
an alternative reference position to compare with the BCG loca- 
tion (see Section 4.1). To minimize the impact of Poisson noise, we 
constructed a smoothed 0.5-2.0 keV image for each cluster using 
the wavelet decomposition method of Vikhlinin et al. (1998). The 
X-ray peak was determined as the highest pixel nearest to the X-ray 
centroid, ignoring any non-cluster point sources, which were iden- 
tified by being detected on only the smallest 1-2 wavelet scales. We 
have neglected the (small) uncertainty in the measurement of the 
X-ray centroid and peak, and therefore also in the corresponding 
projected distance to the BCG, since this is insignificant compared 
to the uncertainties on other parameters. 



2.2 BCG identification and properties 

The positions of the BCG optical peaks were obtained from sev- 
eral sources: near-infrared imaging from the Hale 200-in telescope 
at Palomar Observatory^, described in detail in Stott et al. (2008); 
near-infrared data from the CTIO 4-m telescope^ (Hamilton-Morris 
et al., in preparation); optical imaging from VLT/FORS2'* (PI 
Edge), as well as the Digitized Sky Survey. Photometry of these 
imaging data were supplemented by spectroscopic confirmation of 
membership, using spectra from several sources (see below). The 
BCG was identified as the brightest galaxy in the central few hun- 
dred kpc of the cluster In a few cases the BCG could not be iden- 
tified unambiguously, but these appear to be merging clusters in 



^ The Hale Telescope at Palomar Observatory is owned and operated by 
the California Institute of Technology. 

^ Operated by the Association of Universities for Research in Astronomy 
Inc. 

* Operated by the European Southern Observatory 
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Figure 1. Left: A comparison of the residuals from the Ascasibar & Diego (2008) model, normalized by the measurement errors on each point, as a function 
of scaled radius for both the gas density and temperature data; each line represents a different cluster. Dotted lines mark ±2(7 and the marginal distribution of 
the data is depicted as a kernel smoothed density estimate, compared to a Gaussian of unit variance. Right: A comparison of gas fractions measured within 
both rsoo and r25oo for the four clusters in common with the sample of Vikhlinin et al. (2006); the dotted hne marks the locus of equality. 



which, regardless of the ambiguity, the X-ray/BCG projected offset 
is large. 

Each galaxy was classified according to the presence or ab- 
sence of Hq emission lines, as determined from the literature 
Crawford et al. (1999), Sloan Digital Sky Survey archival spectra 
(Adelman-McCarthy et al., 2008) and VLT F0RS2 spectra for RE- 
FLEX BCGs (PI Edge). For two clusters (Abell 115b and RXC 
J221 1.7-0350) there was no available optical spectrum to determine 
the emission line status of the BCG. Hq emission is a reliable indi- 
cator of the presence of central cold material, second only to direct 
detection of molecular gas at sub-mm wavelengths (Edge 2001). In 
addition, we have classified each BCG according to the presence of 
radio emission detected in the Very Large Array (VLA) FIRST Sur- 
vey (Becker, White & Helfand 1995), NVSS (Condon et al. 1998) 
and/or SUMSS (Bock, Large & Sadler 1999) radio surveys. While 
these three radio surveys have different flux density limits, reso- 
lution and frequency, they nevertheless provide a uniform limit of 
3mJy at 1.4 GHz. 



3 X-RAY CLUSTER MODELLING 

The deprojected gas temperature and density profiles were 
fitted jointly to the phenomenological cluster model of 
Ascasibar & Diego (2008), using the statistic, in order to 
determine the gravitating mass profile and parametrize the gas 
thermodynamic structure (Sanderson, in preparation). The best-fit 
cluster model was also used to estimate rsoo , the radius enclosing 
a mean overdensity of 500 with respect to the critical density 
of the Universe at the cluster redshift. Any spectral bins lying 
within 5 kpc of the centre were excluded from the fit, owing to 
systematic deviations from the model which are expected to be 
significant on small scales (Ascasibar & Diego 2008). For the 
clusters ZwCl 0839.9+2937 and Abell 1758N it was necessary 
to exclude the central bins from the fit (lying at 10 and 30 kpc, 
respectively) due to significant deviations from the best-fit model. 

The left panel of Fig. 1 shows the residuals from the model 
for both the gas temperature and density, as a function of scaled 



radius. It is clear that the model provides a good description of 
the data, with no indication of any significant systematic devia- 
tions. The marginal distribution of the residuals is consistent with 
the Gaussian distribution expected from the measurement errors 
alone, with only a small shift discernable in the case of the temper- 
ature profile- indicating that on average the model tends to slightly 
(~0.3cr) underpredict the data. There is particularly good agree- 
ment between the data and the model in the gas density at all radii, 
including 0.04r5oo, where we use the logarithmic gradient to quan- 
tify the strength of cooling (Section 4.3). It should be noted that the 
gas density profiles extend only to the penultimate spectral bin, ow- 
ing to the non-trivial volume element associated with the outermost 
annulus in the deprojection (Sanderson et al. 2006). 



As an additional test of the model, the right panel of Fig. 1 
plots the derived gas fraction (/gas) within both 7-2500 and rsoo 
compared to measurements made by Vikhlinin et al. (2006), for the 
four clusters common to both samples. We calculate errors on /gas 
and all other derived quantities directly, using the median abso- 
lute deviation of 200 bootstrap resamplings of the input tempera- 
ture and density profile data. There is generally reasonable agree- 
ment in /gas, particularly for the measurements within r2soo, and 
no indication of serious systematic trends, although both Abell 907 
and Abell 2390 are significant outliers, within rsoo. This disagree- 
ment points to the importance of systematics associated with cluster 
modelling, for example differences in the number and location of 
radial bins, as well as the choice of mass profile parametrization. 
Nevertheless, we emphasize the close match between our cluster 
model fits and the data, as seen in residuals from the model. Fur- 
thermore, in Fig. 7 (Section 4.4) we demonstrate good agreement 
between rsoo compared to measurements made by Maughan et al. 
(2008), where we address the issue of bias in X-ray mass estimates 
in relation to cluster dynamical state. 
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4 RESULTS 

4.1 X-ray/BCG offset 

A number of recent studies have identified the projected separa- 
tion between the BCG and the peak of the cluster X-ray emis- 
sion as in important parameter in understanding the activity of 
the BCG, both in terms of star formation (e.g. Edwards et al. 
2007; Bildfell et al. 2008; Cavagnolo et al. 2008) and AON radio 
emission (e.g. Cavagnolo et al. 2008; Mittal et al. 2008). However, 
while the X-ray peak can accurately identify the focus of a cool 
core, this location may not always lie at the 'centre' of the hot gas, 
as defined for the purposes of conducting a radial profile analy- 
sis. This is particularly true for non-CC clusters, which have much 
broader and flatter X-ray cores. We therefore investigate the pro- 
jected offset between the BCG and both the peak and centroid of 
the X-ray emission. We assess the relative merits of both definitions 
in the context of quantifying cluster dynamical state in Section 4.3, 
but turn our attention initially to the projected offset in terms of the 
centroid. 

Fig. 2 shows the distribution of X-ray centroid/BCG projected 
offsets for the sample, in units of arseconds, kiloparsecs and rsoo- 
The distributions are roughly lognormal, with medians of 13kpc 
/ 1.0 per cent of rsoo. Also plotted are the separate distributions 
for BCGs with and without detected Ha line emission, and a clear 
difference is apparent: the BCGs with line emission all have small 
projected offsets from the X-ray cluster centroid, within a maxi- 
mum of 15 kpc / 1.5 per cent of rsoo. Conversely, BCGs with no 
detected Hq emission lines have preferentially larger X-ray cen- 
troid/BCG projected offsets (median values = 35 kpc / 2.8 per cent 
of rsoo)- This confirms that close proximity of the BCG to the X-ray 
cluster centroid is a pre-requisite for BCG star formation, in agree- 
ment with recent work (Edwards et al. 2007; Bildfell et al. 2008; 
Rafferty et al. 2008). 

For comparison. Fig. 3 shows the distribution of projected off- 
sets for BCGs with and without detected radio emission. As with 
the Ha emitters in Fig. 2, the clusters with radio detected BCGs 
have smaller X-ray/BCG offsets. However, a separate peak of radio 
detected BCGs with intermediate projected offsets is also present. 
For the sample as a whole, the interquartile range of X-ray/BCG 
projected offsets is 0.0035-0.04rsoo (see bottom right panel of 
Fig. 5), demonstrating that BCGs in X-ray selected clusters are gen- 
erally located close to the cluster X-ray centroid; all but two (97^4 
per cent) lie within 0. 15r5oo. 
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Figure 2. Kernel-smoothed probability density estimates of the distribution 
of projected offsets between the X-ray cluster centroid and the BCG. The 
raw values are shown at randomly 'jittered' y-axis positions. The distribu- 
tions for BCGs with and without detected Ha emission are also indicated. 



Emission lines No lines 

Radio emission 21(0.33) 10(0.16) 

No radio 2 (0.03) 30 (0.48) 



Sum 

31 (0.51) 

32 (0.49) 



Sum 



23 (0.36) 40 (0.64) 63 (1) 



Table 1. A contingency table of the BCG distribution according to their Ha 
and radio emission status, with fractions of the total (63) given in brackets. 
Note that the two clusters with unknown emission line status are excluded. 



4.2 BCG Ha and radio emission 

Table 1 sunmiarizes the Ha and radio emission status of the 63 
clusters for which both classifications are available. Roughly half 
the galaxies are passive (30/63), with no detected Ha or radio emis- 
sion. The radio active fraction is 49±7 per cent, which is larger 
than the value of > 30 per cent found in clusters selected opti- 
cally (Best et al. 2007) and in X-rays (Lin & Mohr 2007), possi- 
bly due to the high average mass of the clusters in the LoCuSS 
sample. Similarly, the fraction of Ha emitting BCGs (37 ± 7 per 
cent) is somewhat higher than the value of 27 percent for the X-ray 
selected cluster sample of Crawford et al. (1999), which is itself 
roughly double the fraction in optically-selected clusters (13 per 
cent; Edwards et al. 2007). A Pearson's test strongly rejects the 
null hypothesis that the presence or absence of detected Ha emis- 
sion is independent from the presence or absence of radio emission, 
with a value of 25.7 for 1 degree of freedom and a correspond- 



ing p-value of 4.0 x 10~^. This demonstrates the close coiTespon- 
dence between BCG activity in terms of star formation and black 
hole accretion. 

4.3 Cluster cool core status 

In order to connect the properties of the BCG to those of its host 
cluster, we focus on the cooling state of the intracluster medium 
(ICM). A key signature of cool core clusters is a cuspy gas density- 
and hence surface brightness- profile, which can be more accu- 
rately measured than the associated decline in temperature. Fur- 
thermore, the gradient of the gas density profile on small scales 
progressively steepens with increased cooling, even before a signif- 
icant cool core is established (Ettori & Brighenti 2008). Therefore, 
following Vikhlinin et al. (2007), we use the parameter a, defined 
as the logarithmic slope of the gas density profile at 0.04r5oo, to 
quantify the extent of cluster gas cooling on the approximate scale 
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Figure 3. Kernel-smoothed probability density estimates of the distribution 
of projected X-ray /BCG offsets, normalized to rsoo ■ The data are separated 
according to whether or not radio emission is detected from the BCG and 
the raw values are shown at randomly 'jittered' y-axis positions. 
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Figure 4. a (see Section 4.3) versus central gas entropy, from the best- 
fit cluster model. The vertical dashed line marks 30 keV cm^ . Point styles 
indicate the presence or absence of Ho emission lines in the BCG and the 
colour shows whether significant radio emission is detected from the BCG. 
The zero value for Abell 1423 is plotted as an arrow from its upper bound. 



of any central galaxy, with more negative values implying stronger 
cooling. 

We calculate the logarithmic slope, a, from the best-fit 
Ascasibar & Diego (2008) cluster model and plot it, in Fig. 4, 
against the central gas entropy, which is a direct indicator of the 
cooling state of the ICM. However, this measurement involves ex- 
trapolation of the model to zero radius, unlike a, which is evaluated 
at a directly accessible radius, which also lies in a region where the 
model provides a good fit to the density profile (Fig. 1, left panel). 
The strong correlation in Fig. 4 demonstrates the effectiveness of a 
as a simple means of quantifying cool-core strength. 

It can be seen from Fig. 4 that there is a gap in the en- 
tropy values which separates the clusters into two regimes, with 
all the Hq emitting BCGs lying in clusters with central entropy 
<30keVcm^, as discovered by Rafferty et al. (2008) and subse- 
quently confirmed by Cavagnolo et al. (2008). Furthermore, there 



are a number of clusters with a central gas entropy of ~ 1 ke V cm^ , 
with the remainder grouped around '^l 00-200 keVcm^ as also 
found by Cavagnolo et al.. This threshold entropy level can be un- 
derstood as being determined by the physics of thermal conduction 
(Voit et al. 2008), which can also account for the bifurcation into 
CC and non-CC populations (e.g. Guo, Oh & Ruszkowski 2008; 
Sanderson et al. 2009). 

Given the greater ease and reliability with which a can in gen- 
eral be measured, we use this as the parameter of choice for quan- 
tifying cool core strength. We now turn to the issue of how cluster 
cooling is related to the X-ray/BCG offset. Fig. 5 shows the varia- 
tion of a with X-ray/BCG projected offset (scaled by rsoo) defined 
using both the X-ray centroid (left panel) and peak (right panel). A 
clear trend is evident, with the strongest cooling cores found only 
in clusters with small projected X-ray/BCG offsets. Vikhlinin et al. 
(2007) suggest a < —0.7 for strong cooling flows, which lies close 
to the gap in the distribution of values. However, it is clear that only 
those clusters with a < —0.85 and a projected offset < 0.02r5oo 
have detected Hq line emission (as depicted by the shaded box), 
which is indicative of genuine gas condensation and subsequent 
star formation associated with unchecked cooling. Of the 24 sys- 
tems in this shaded region, only 1 is passive, and 21 have both Ha 
and radio emission; there are no Ha emitting BCGs outside the 
shaded region. 

While the two panels in Fig. 5 show similar results, the cor- 
relation between a and the X-ray centroidfBCG projected offset 
is stronger and shows fewer outliers: the Kendall rank correlation 
coefficient gives r = 0.48 (p-value = 1.3 x 10^*) for the offset 
with respect to the centroid and r = 0.42 (p-value = 6.7 x 10^^) 
for the peak. Nevertheless, a number of prominent outliers are 
still present in the left panel of Fig. 5 and are labelled in both 
panels. In particular Abell S0592 and RXC J0232.2-4420, which 
have relatively large X-ray/BCG offsets, despite hosting strong 
cool cores. Abell S0592 has an obvious subclump (excluded from 
the X-ray analysis) which is associated with the BCG; the main 
peak X-ray centroid is located close (<10kpc) to another galaxy. 
RXC J0232. 2-4420, on the other hand, has a relatively regular and 
single-peaked morphology, but with 2 bright galaxies of which the 
less luminous is located close to (<10kpc) the X-ray centroid. 

The only passive cluster with a steep a and small projected 
X-ray/BCG offset is Abell 1423. While there is a extended radio 
source very close to the BCG, inspection of the FIRST Survey im- 
age indicates that this is an unrelated narrow angle tail radio galaxy 
seen in projection and that the BCG is undetected to the FIRST 
depth. No optical lines were detected by Crawford et al. (1999) but 
the spectrum was not of the highest quality so more detailed follow- 
up of this BCG would be of great interest. The two remaining la- 
belled points are the clusters RXC J0220.9-3829 and Abell 3088, 
which are the only other systems in the shaded box with BCGs that 
have no significant radio emission. 

While it is clear that Ha emitting BCGs are all concentrated in 
the bottom left corner of Fig. 5 (particularly the left panel), the dis- 
tribution of passive BCGs (blue circles, with no Ha or radio emis- 
sion) follows a continuous strong trend: a test with the Kendall rank 
correlation coefficient for the centroid offset case gives a value of 
T 0.53 (p- value = 2.0x 10^'') in favour of the null hypothesis that 
the two quantities are uncorrelated. The equivalent Kendall coef- 
ficient for the peak offset case is r 0.48 (p-value = 1.0x10"^). 
This suggests a progressive weakening of central gas cooling with 
increasing cluster disruption, as measured by the X-ray/BCG pro- 
jected offset, even for clusters with no evidence of condensation 
from the hot phase. 
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Figure 5. a (see Section 4.3) plotted against the projected offset between the BCG and the X-ray centroid (left panel) and the X-ray peak (right panel). The 
shaded box highlights the region occupied by line emitting galaxies in the left panel. 
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Figure 6. Cluster gas fraction within rsoo as a function of projected offset between the BCG and the X-ray centroid (left panel) and the X-ray peak (right 
panel). The horizontal shaded region depicts the Universal baryon fraction as measured from the five year WMAP data (Dunkley et al. 2009). Outliers identified 
in Fig. 7 ai'e labelled: see text for details. 



4.4 Cluster gas fraction and dynamical state 

To explore the connection between X-ray /BCG offset and cluster 
dynamical state, we plot in Fig. 6 the gas fraction (/gas) within rsoo 
as a function of projected X-ray/BCG offset, with the Universal 
baryon fraction (/b) from the five year WMAP data (Dunkley et al. 
2009) also indicated. As in Fig. 5, the relation is shown for the off- 
set in terms of both the X-ray centroid (left panel) and peak (right 
panel). A test with the Kendall coefficient for the centroid offset 
case gives t = 0.31 (p- value = 3.0 x 10^^) in favour of the null hy- 
pothesis that /gas and the X-ray/BCG offset are uncorrelated. The 
equivalent Kendall coefficient for the peak offset case is r = 0.24 
(p-value = 5.0x10""'). In contrast, there is no significant correla- 
tion between the X-ray/BCG offset and the cluster mean tempera- 
ture or total mass, which might otherwise account for such a trend 
with /gas (cf. Sanderson et al. 2003; Vikhlinin et al. 2006) 

All of the five clusters statistically consistent with having gas 
fraction values in excess of the WMAP5 /b upper bound are pas- 
sive systems and all have projected X-ray/BCG offsets exceed- 
ing 0.03r5oo- all are also identified as outliers in Fig. 7 (see be- 



low). The cluster gas fraction within rsoo is typically ~90 per cent 
of the cosmic mean (e.g. Grain et al. 2007) and the baryon frac- 
tion in massive clusters is roughly 10 per cent higher than /gas 
(e.g. Lin, Mohr & Stanford 2003; Gonzalez, Zaritsky & Zabludoff 
2007). Therefore such large gas fractions are unrealistic and may 
indicate a failure of the assumption of hydrostatic equilibrium 
in the cluster modelling, as might result from significant dynam- 
ical disturbance (Nagai, Vikhlinin & Kravtsov 2007). These ob- 
jects are therefore likely to be outliers and significant sources of 
scatter in scaling relations. Taking the nominal X-ray/BCG offset 
value of 0.02r5oo used to identify the strongest cool core clus- 
ters (see Fig. 5), the mean gas fraction of the 'small offset' clus- 
ters (^0.02r5oo) is 0.106 ± 0.005 (a = 0.03) compared to 
0.145±0.009 (a = 0.04) for the clusters with an offset > 0.02r5oo. 

The only emission line BCG cluster with a high gas frac- 
tion in Fig. 6 is Abell 2390 (/gas = 0.17 ± 0.01, lying within 
the WMAP5 /b band), which is a hot cluster (9.8 keV) that had 
the highest /gas in the 13 cluster sample of Vikhlinin et al. (2006) 
and which is a known merging system, despite hosting a cool core 
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(Allen, Ettori & Fabian 2001). As pointed out above, the compar- 
ison between our /gas value and that of Vikhlinin et al. (Fig. 1 , 
right panel) reveals a clear discrepancy within rsoo. Neverthe- 
less, the agreement is better within 7-2500, and our rsoo value of 
1437±68 kpc is very close to the value of 1416±48 kpc calculated 
by Vikhlinin et al.. Furthermore, it can be seen from inspection of 
figure 12 in Vikhlinin et al. (2006) that their gas fraction profile for 
Abell 2390 rises sharply just outside rsoo, reaching values consis- 
tent with our measurement, albeit at a slightly larger radius. 

To explore the effectiveness of the projected offset between 
the BCG and the X-ray centroid in quantifying the cluster dy- 
namical state, we turn to a comparison between rsoo determined 
from two different methods. Fig. 7 shows the comparison be- 
tween rsoo as calculated from the Ascasibar & Diego cluster model 
and that measured by Maughan et al. (2008) using a scaling rela- 
tion based on the X-ray equivalent to the Compton y-parameter 
(Ix; Kravtsov, Vikhlinin & Nagai 2006). Yx is simply the total 
thermal energy of the hot gas within a given radius, and acts 
as a robust mass proxy, even for dynamically disturbed clusters 
(Kravtsov et al. 2006; Poole et al. 2007; Maughan 2007). The 30 
clusters common to both samples are split into four subsamples de- 
fined by quartiles of the projected X-ray centroid/BCG offset. In 
general, the agreement is excellent, particular for the clusters in the 
50th percentile of X-ray centroid/BCG projected offset (top panels 
of Fig. 7), however six prominent outliers are visible and have been 
labelled. 

The Ascasibar & Diego cluster model assumes hydrostatic 
equilibrium (HSE) to derive rsoo, whereas rsoo inferred from Yx 
is only sensitive to the total thermal energy of the gas. A compari- 
son between the two values therefore reveals the extent to which the 
assumption of HSE applies for any given cluster. For five of the out- 
liers in Fig. 7 the Ascasibar & Diego underpredicts rsoo compared 
to the Yx-based measurement of Maughan et al., consistent with 
the tendency for HSE-based measurements to underestimate the 
cluster mass (Nagai et al. 2007). These outliers all fall in the high- 
est two quartiles of X-ray centroid/BCG projected offset, further 
supporting the view that this offset acts as an indicator of dynam- 
ical disturbance (cf. Katayama et al. 2003). The remaining cluster, 
Abell 267, has a complex mass distribution, based on the gravita- 
tional lensing analysis of Smith et al. (2005), which has a highly 
elliptical morphology. This could result in the HSE-based X-ray 
analysis overestimating the total mass, as suggested by the compar- 
ison of rsoo in Fig. 7, as well as its low /gas (0.041 ±0.018; Fig. 6), 
given its large mass (Msoo=1.85 ± 0.70 X lO^^M©). 



5 SUMMARY AND DISCUSSION 

We have shown that the projected X-ray/BCG offset is highly cor- 
related with the strength of cooling in the host cluster core, as mea- 
sured by the logarithmic slope of the gas density profile at 0.04rsoo 
(a; Section 4.3) and that the gas fraction is systematically larger in 
clusters with large offsets. In particular, the use of the X-ray cen- 
troid yields a stronger correlation in both these relations, compared 
to using the X-ray peak to calculate this offset. 

Our results also demonstrate that the activity of BCGs is 
closely related to the properties of their host cluster, and specifi- 
cally the proximity to and strength of any cool core. The behaviour 
of the Hof emitting galaxies is especially striking in that all 23 of 
them in our sample of 65 BCGs are found close to (< 0.02rsoo) a 
strong (a < —0.85; Section 4.3) cool core, with only 1 other non- 
emission line galaxy meeting these same criteria. This is consistent 
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Figure 7. The comparison between rsoo calculated from the 
Ascasibar & Diego (2008) cluster model and that measured by 
Maughan et al. (2008) for the 30 clusters common to both samples. 
The clusters have been split into quartiles according to their projected 
X-ray centroid/BCG offset, increasing from the top left to bottom right, 
as indicated by the range in the strip for each panel (0.00187-0.00467; 
0.00467-0.0164; 0.0164-0.0565; 0.0565-0.34 rgoo). The dashed lines 
indicate the locus of equality. 

with the work of Rafferty et al. (2008) and Cavagnolo et al. (2008) 
who find that star formation only occurs once the hot gas entropy 
falls below a critical threshold of ~30keVcm^. Voit et al. (2008) 
show that such a threshold can be understood if energy output from 
AGN is coupled to the cooling gas via thermal conduction, a pro- 
cess which is also capable of thermally stabilizing no«-cool core 
clusters (Sanderson et al. 2009). 

Notwithstanding the effects of conduction or galaxy feedback, 
it is clear that star formation is ultimately able to take place when 
the BCG coincides closely with a strong cool core. As a conse- 
quence, additional enrichment of the ICM is possible in the vicin- 
ity of the BCG, which may explain the centrally peaked metal- 
licity profiles of cool core clusters (Leccardi & Molendi 2008; 
Sanderson et al. 2009) and the fact that the lowest entropy gas is 
also the most enriched (Sanderson et al. 2009). The influence of 
the BCG on the ICM may also account for the increased varia- 
tion in gas properties (e.g. metallicity and entropy) on small scales 
(~0.02r5oo; Sanderson et al. 2009). 

The strong connection between line emission and radio emis- 
sion (Section 4.2) suggests that star formation and AGN accretion 
are fuelled from the same source. However, a significant fraction 
(23 per cent; 7/31) of the radio emitting BCGs occupy a nan'ow 
range of relatively large offsets 0.032-0.079r5oo and have no Ha 
emission, suggesting a different origin. These may be cases where 
recent merger disruption has triggered a brief burst of AGN activity. 

If the X-ray/BCG projected offset serves to measure the dy- 
namical state of the cluster, with more disturbed systems having 
larger values, then its strong correlation with the steepness of the 
gas density profile (a) implies that cool core strength progressively 
diminishes in more dynamically disrupted clusters. Such a trend 
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would be expected if cluster mergers were capable of erasing cool 
cores. There is a notable dearth in Fig. 5 of clusters with relatively 
flat density profiles and small projected X-ray/BCG offsets. Such 
cases ought to exist if AGN outbursts alone were capable of reheat- 
ing cool cores, or if strong preheating prevented their formation 
altogether (McCarthy et al. 2008). On the other hand, if such pro- 
cesses had a disruptive impact on the core X-ray morphology, then 
this might shift the centroid as measured on ~arcmin scales enough 
to explain the trend. 
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